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We report on preliminary measurements of time-dependent CP asymmetries in neutral B decays to 
CP eigenstates with transition amplitudes that are dominated by penguin-type loops. The results 
are obtained from a data sample of up to 227 million T(45) — > BB decays collected with the BABAR. 
detector at the PEP-II asymmetric-energy B-meson Factory at SLAC. The amplitudes of the effective 
mixing-induced CP asymmetries, sin2/3 c ff , are derived from decay-time distributions of events in which 
one neutral B meson is fully reconstructed in a final state without charm and the other B meson is 
determined to be either a B° or B° from its decay products. 



1 Introduction 

Since the precise measurement of sin2/3 from 
the time-dependent CP- asymmetry analyses 
of b — > ccs decayflEl considerable effort 
is spent at the B-Factory experiments to 
determine the corresponding quantity from 
penguin-loop-dominated (charmless) B de- 
cays ("s-penguin decays"), such as B° — > 
(/)K , r/'K , f {980)K°, S-w&ve dominated 
K + K~K° and tt°K°. Due to the large vir- 
tual mass scales occurring in the penguin 
loops, additional diagrams with heavy par- 
ticles in the loops and new CP-violating 
phases may contribute to the decay ampli- 
tudes. The measurement of CP violation 
in these channels and the comparison with 
the charmonium reference value is therefore 
a sensitive probe for physics beyond the 
Standard Model (SM). Due to the different 
non-perturbative strong-interaction proper- 
ties of the various s-penguin decays, the NP- 
induced deviation from sin2/3 is expected to 
be channel-dependent. 

There is a variety of reasons why spe- 
cific contributions of New Physics (NP) are 
expected at or below the TeV scale. Among 
these are the gauge hierarchy problem of the 
SM Higgs sector, Baryogenesis, Grand Unifi- 
cation of the gauge couplings and the strong 



CP problem. If the NP is not flavor blind 
it should be present in the flavor se ctor , and 
detectable in loop-induced B decays.^] How- 
ever, the present insignificance of NP effects 
already sets strong boundary conditions on 
the generic flavor structure of the NP, hence 
favoring minimal flavor- violating schemes.^ 

Assuming penguin dominance and ne- 
glecting CKM-suppressed amplitudes, s- 
penguin decays carry approximately the same 
weak phase as the decay B° — > J/ipK®. Their 
mixing-induced CP- violating parameters are 
expected to b d 1 1 —ijf x sin2/3 ~ —r]f x 0.7 
in the SM a , where rjt is the CP eigenvalue of 
the final state /. 

This simple relation is altered by the 
presence of a suppressed V us V* b penguin di- 
agram in all s-penguin decays, and a corre- 
sponding color-suppressed tree-level diagram 
in three-body decays (such as K + K~ K°) or 
resonant decays with non-ss flavor contri- 
butions (like n'K°, foK n and tt K°). As 
a consequence, only an effective sin2/3 c ff = 
—r\f x S is determined. Naive dimensional 
arguments may allow us to derive coarse 
estimates for the mode-specific theoretical 
systematics on the amplitude level. We 



"With these assumptions, the expected phase shift is 
2/3 cH - 2/3 = -axg[(V t6 V t * V* b V cs ) 2 \ ~ 2.1°. 
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find an uncertainty of 0(5%) (correspond- 
ing to a 2/3 e ff shift of about 2°) for the 
golden mode (f>K°, 0(10%) for the silver- 
plated modes rj'K , f K° or K + K~K°, and 
0(20%) for the bronze-plated modes that 
have no ss component (Tr a K° and similar). 
Computations based on QCD factorization 
generally find smaller deviations.^ Bounds 
on 2/3 c ff — 2/3 can be derived from SU(3) 
flavor-symmetry relations^ which are how- 
ever weak at present. 

We present in this summary preliminary 
measurements of CP-violating asymmetries 
in all of the above mentioned s-penguin de- 
cays. The data used for these measure- 
ments were accumulated with the BABAR 
detected at the PEP-II asymmetric-energy 
e + e~ storage ring at SLAC. The sample con- 
sists of 227 x 10 6 BB pairs (209 x 10 6 for 
70(980)^°) collected at the T(AS) resonance. 
Detailed descriptions of the individual anal- 
ysespresented here are given in RefsPl. In 
Rcf. * we describe the silicon vertex tracker 
(SVT) and drift chamber used for track and 
vertex reconstruction, the Cherenkov detec- 
tor (DIRC), the electromagnetic calorime- 
ter (EMC), and the instrumented flux return 
(IFR). 

With At = tcp ~ ttag defined as the 
proper time interval between the decay of the 
fully reconstructed B^ P and that of the other 
meson B^ g , the time-dependent decay rate is 
given by 



-|At|/T 

At 



l + QtagS , sin(Am d At) 



- Q tag Ccos(Am d At) , (1) 



/Q tag (Ai) 



where Qtag = 1( — 1) when the tagging me- 
son B t ° ag is a B°(B°), r is the mean B° life- 
time, and Arrtd is the B°B° oscillation fre- 
quency. The parameter 5* is non-zero if there 
is mixing-induced CP violation, while a non- 
zero value for C would indicate direct CP vi- 
olation. At least two interfering decay ampli- 
tudes with different weak and strong phases 



are required to obtain C ^ 0. Hence we do 
not expect to observe significant direct CP 
violation in any of the channels studied if the 
assumption of penguin dominance is correct. 



2 Analysis Techniques 

We reconstruct signal candidates from combi- 
nations of a K® candidate decaying to tt + it~ , 
and a <fi — > K + K~ , rf — ► (r;7r + 7r~ , p°j), 
Jo (980) — > 7r + 7r - , K + K~ (excluding </>), or a 
7T° candidate. We also reconstruct the mode 
B° — > 4>K®, where a K Q L candidate is identi- 
fied either as an unassociated cluster of en- 
ergy in the EMC or as a cluster of hits in 
the IFR. For the B° -> r)'K Q g decay, also 
K° — > 7r°7r° and 7/ — > (77, 7r + 7r~7r ) decays 
are considered so that the reconstructed sig- 
nal amounts to 46% (not containing the de- 
tector acceptance) of the total r]'K® produc- 
tion. 

We use strong particle identification from 
the DIRC, the tracking system and the EMC, 
to identify pions and kaons in the final state. 
Two kinematic variables are used to dis- 
criminate between signal- 5 decays and com- 
binatorial background: the difference AE 
between the center-of-mass (CM) energy of 
the B candidate and the CM beam en- 
ergy, and the beam-energy-substituted B- 
candidate mass wies, defined in the labora- 
tory frame. In the tt°K® analysis, these vari- 
ables are replaced by the unconstrained B- 
candidate mass m rec = \qb\-, and the missing 
mass m m i ss = \q e + e - — Qb\, where q e + e - is the 
four-momentum of the e + e~ system and qs 
is the four-momentum of the B candidate af- 
ter applying a B°-mass constraint. The miss- 
ing mass has a slightly better resolution than 
Trass and, by construction, m rcc and TO m ; ss 
have a vanishing linear correlation coefficient. 

Continuum e + e~ — * qq (q = u,d,s,c) 
events are the dominant background. To 
enhance discrimination between signal and 
continuum, we use a multivariate analyzer 
(MVA), which is a Fisher discriminant or a 
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neural network (NN), to combine four dis- 
criminating variables: the cosine of the an- 
gle between the CM B direction and the z 
axis (along the beam direction), the cosine of 
the angle between the thrust axis of the B 
candidate and the z axis, a nd th e topological 
monomials Lq, L2 (see RefsPIl). 

Monte Carlo (MC)-simulated events are 
used to study the background from other B 
decays. The charmless modes are grouped 
into different classes with similar kinematic 
and topological properties. In the K + K~K° 
analysis exclusive b — > c decays with the same 
final state as the signal are vetoed by cuts on 
their invariant masses. 

The time difference At is obtained from 
the measured distance between the z posi- 
tions of the B^p and B® decay vertices, and 
the boost /?7 = 0.56 of the e + e~ system. 
Due to the K® lifetime, the At determina- 
tion is particularly challenging for the Tr°K® 
final state. To select valid candidates for this 
mode, at least 4 SVT hits are required for 
each of the decay pions (about 60% of the 
events) . While the At information of the re- 
jected events is not used, they contribute to 
the measurement of the C coefficient. We 
determine At from a global constrained fit to 
the entire T(45) B°B° decay tree that 
takes the information on the beam energy 
and the position of the interaction point into 
account. To further improve the At resolu- 
tion, the sum of the two measured B° life- 
times is constrained to 2t with an uncertainty 
of \P1t in the fit. 

A new NN-based algorithnJEl 

has been 

developed to determine the flavor of the B® ag . 
In addition to the information exploited by 
the algorithm described in RefP, it uses low- 
momentum electrons, A — > pn decays, and 
additional correlations among identified kaon 
candidates. The resulting effective tagging 
efficiency amounts to Q ~ 30.5%. The tag- 
ging efficiencies, mistag probabilities, mistag 
biases and also the At resolution parame- 
ters are determined from a fit to fully recon- 



structed B decays to charm. 

We use unbinned extended maximum- 
likelihood fits to extract the event yields and 
the CP parameters defined in Eq. QJ. In 
these fits, the likelihood for a given event is 
the sum of the signal, continuum and the B- 
background likelihoods, weighed by their re- 
spective event yields. We use as fit variables 
At, AE, m ES , the MVA (all modes) as well 
as the resonance mass and decay angle (<fiK® 
and / (980)JC°). In the case of <pK° the B- 
mass kinematic constraint is necessary to de- 
termine the K® momentum so that wes can- 
not be exploited in the fit. The signal refer- 
ence shapes for the fit variables are obtained 
from MC simulation, validated with control 
samples from fully reconstructed B decays to 
charm of similar topology. As many back- 
ground shape parameters as possible are si- 
multaneously determined by the fit to reduce 
the use of prior assumptions. 

The K + K~K® final state has unknown 
a priori CP content so that the CP-even frac- 
tion, /oven, has to be extracted experimen- 
tally. Two different methods are applied, of 
which the first one (described below) pro- 
vides the main result, and the second is a 
cross check. We study the K + K~ helicity 
angle distribution \A B a^ K + K - K o s (cos#)| 2 = 
Sf=o,i.2(^) x Pe(cos0), where P t are Leg- 
endre polynomials of order I. The aver- 
age moments (Pi) are computed by means 

of background-subt ract ed signal weights re- 
in J)| 

turned by the fitM^ 1 We find / eTcn = 

A s/( A s + a p) = °- 89 ± °- 08 ± °- 06 ' where 
the first error is statistical and the sec- 
ond systematic, and where the S- and P- 
wave amplitude-squared are given by A 2 S = 
V2(Po) - V / 572(- p 2) and A% = <ftj2(Fi). 
Assuming isospin invariance, fe vpn can also 
be obtained from the raticQSl 2T(B + — > 
K + K°K°)/T(B° -> K+K~K°), for which 
we find 0.75 ± 0.11 (statistical error only), in 
agreement with the result from the moments 
analysis. 
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Mode Yield sin2/? eff C 



B - 




(<f>K° s ) 

mi) 


114 ± 12 
98 ± 18 


+0.50 ± 0.25 ig;g| 


0.00 ±0.23 ±0.05 


B - 


- V'K° S 




819 ±38 


+0.27 ±0.14 ±0.03 


-0.21 ±0.10 ±0.03 


B - 


■* / (980)tf° 




152 ± 19 


+0.95^;^±0.10 


-0.24 ±0.31 ±0.15 


B - 


■+ K+K-K° s 




452 ± 28 


+0.55 ± 0.22 ± 0.04 ± 0.11 C p 


+0.10+0.14+0.06 


B - 






300 ± 23 


+0.35 ±^ ±0.04 


+0.06 + 0.18 + 0.06 



Table 1. Preliminary fit results for the different s-penguin modes. The second column gives the observed 
signal yield, the third the effective sin2/3 c fj = — r/f X S value, and the last column the result for the C parameter. 
The first errors given are statistical and the second systematic. The third error given for sin2/3 c ff (K+ K~ Kg) 
is due to the uncertainty on the CP-even fraction. 



3 Results 

The fit results for the event yields and CP pa- 
rameters are given in Tabled I n the case of 
K + K~ K®, sin2/3 e ff and / evcll (approximately 
related through sin2/3 c g ~ S/(2f even — 1)) 
are simultaneously determined by the fit. 
The individual sin2/3 c ff results for the <fiK® 
(0.29 + 0.31) and <j>K° (1.05 + 0.51) subsam- 
ples agree within statistical errors. We have 
also measured the charge asymmetry in the 
decay B + — > <fiK + for which we find the null 
result A^ K+ = (r^if") - T(4>K+))/£ = 
+0.054+ 0.056+ 0.012. 

All measurements reported here are 
statistics limited. The systematic errors 
are dominated by the opposite CP back- 
ground and the reference shape modeling 
((t)K°), the fit bias (rfK°, ,f (980)Jf° and 
K + K~ K^), the interference with other res- 
onances (/o(980)ifg), and the background 
tagging asymmetries and SVT alignment 
(Tr°Kg). All measurements have in common 
the uncertainty due to the unknown CP vio- 
lation on the tag-side. 

The individual sin2/3 ff results are in 
agreement with the charmonium value,^^ 
with the exception of rfK®, which exhibits 
a 3. Oct discrepancy (neglecting theoretically 
motivated deviations). Overall, the results 
are in agreement with but more precise than 
the correspo ndin g numbers from the Belle 
CollaboratiorPSI. 
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Figure 1. Compilation of preliminary results on 
Sin2/3f e ff) from charmonium and s-penguin decays. 
The averages, indicated by the shaded bands, are ob- 
tained ignoring theoretical uncertainties. 



A graphical compilation of the results is 
given in Fig. ^ We find the s-penguin av- 
erage sin2/3 eff = 0.42 ± 0.10 (C.L. = 0.40), 
where theoretical uncertainties are ignored. 
The difference with the charmonium refer- 
ence amounts to 2.7 standard deviations (a). 
Adding the aforementioned individual theo- 
retical uncert ainti es, and treating them as 
allowed ranges pi the significance of the dis- 
crepancy drops to 2.0cr. 

None of the modes studied exhibits non- 
zero direct CP violation. The s-penguin av- 
erage is C = -0.081 ± 0.073 (C.L. = 0.43). 
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4 Summary 

We have presented preliminary results on 
mixing-induced and direct CP violation for 
penguin-loop-dominated B decays, using the 
full Summer 2004 data sample. We have 
observeoP^ and updated the new silver- 
plated mode B° -> / (980)X°. Sophisticated 
vertexinjl^l allows us to measure the mixing- 
induced CP asymmetry in B° — > ^K®, using 
the K® decay into ir + ir~. With the excep- 
tion of B° —> r]'K®, all individual s-penguin 
modes are in agreement with the charmo- 
nium result, but most of them tend to lower 
sin2/3 e ff values so that their average appears 
to be somewhat low. 

With the significant decrease in the sta- 
tistical errors, the B-meson Factories enter 
the era of precision CP-violation measure- 
ments in loop processes. 
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